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1 A novel Lie-group analysis of the (unclosed) mean momentum equation (MME) for turbulent channel and pipe flows 
. yields an analytic multi-layer formula for the entire mixing length profile (hence the mean velocity profile). Simple 
Lie-group arguments lead to three kinds of local dilation-invariant expressions of the mixing length, which we believe 
\ to be the most basic solutions of MME describing each statistically autonomous layer. The first kind - a power-law 
solution - is valid for viscous sublayer, buffer layer, log-layer, and a newly identified central core. Remarkably, the 
second - a defect power-law of form 1 — r rn - describes a recently founded new bulk zone of channel and pipe flows. 
In addition, a simple ansatz describing a relation between invariants of the mixing length and its gradient, is shown 
to describe scaling transition between two adjacent layers. The theory thus proposes a rigorous quantification of the 
multi-layer structure, and identifies three kinds of physically meaningful parameters: group-invariant scaling, layer 
__■ thickness and transition sharpness. Our results are validated by direct numerical simulation (DNS) of channel and 
pipe flows, while extensive contrast with numerous empirical data and the determination of all layers' parameters will 
^> be discussed elsewhere. 

V 

I. INTRODUCTION 

<Z3 

Wall-bounded turbulent flow has received considerable recent attention because of its relatively simple geometry 
and its relevance to a large variety of engineering applications [l| • It enables direct insight into fundamental turbulence 
physics, as well as direct verification by local (such as profiles of mean velocity, Reynolds stress, fluctuation intensities, 
etc.) and integral measurements (such as skin friction and heat transfer). The fundamental challenge is to predict 
the mean flow properties, including the mean velocity profile (MVP), the Reynolds stress, the kinetic energy, etc; 
however, a deductive theory of this kind is still missing. 

A recent model by L'vov, Proccacia and Rodenko (cited hereinafter as LPR) Q is noteworthy, as they deal with the 

■ MME and the balance equations for the kinetic energy and the Reynolds stress. Several physical arguments regarding 
I | energy cascade, dissipation and balance in the average Reynolds stress equation, enable LPR to prescribe the mixing 

^C) ■ length profile which involves two artificial functions inspired from the limited DNS data - one connecting the log-layer 
\Q ' to the wall (e.g. a wall function) and the other to the center (e.g. a bulk function). This two-layer model achieves 

■ a good accuracy in describing the data, but the asymptotic behaviors for the wall and bulk functions are not correct 
Q, and their extensions to describe other flows and other profiles (energy dissipation, kinetic energy, etc) are not 

t-H . obvious. 

t— I ■ Recently, a SED theory [4| has been proposed to reconsider this problem of determining the mean velocity profile 
of canonical wall-bounded flows from first principles. It is known that the mean momentum equation is unclosed, 
and the determination of the mean strain requires a modeling of the Reynolds stress. Prandtl's Q mixing length 
model provides a way to model the Reynolds stress. In the past, simple models for the mixing length were adopted; 
the most famous is the Karman's linear law of the distance to the wall, leading to a logarithmic profile for the mean 
velocity. The log-law has had a tremendous success in predicting the Reynolds number dependence of the friction 
coefficient and in calibrating the computational fluid dynamics (CFD) models of turbulent flows such as K-e or K — ui 
[lj. Nevertheless, the validity of the logarithmic profile, even for simple wall-bounded flows, has been questioned 043; 
owing to the ad hoc nature of the arguments behind its derivation. In Q , we have extended the classical mixing length 
theory of Prandtl and von Karman by proposing a multi-layer formula for the entire profile of the mixing length £^f 
- hence the entire mean velocity profile - written as 
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where y + is the coordinate normal to the wall, r is the radial distance of the pipe, y^ ub , y bu f and r core are parameters 
specifying the thicknesses of the layers (Z c = (1 + r 2 cor ^) x l A is a normalization factor), pi and P2 are parameters 
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describing the sharpness of the transitions between the layers, m is a parameter specifying a geometry-dependent 
scaling for the bulk flow (empirical study shows that m = 4 for channel and m = 5 for pipe [4|), and the superscript 
'+' denotes wall unites normalization by the wall friction velocity u T and viscosity. This formula (TTJ), seemingly 
complex, is in fact composed of a series of simple functions, namely 

e + M ^ (p/ylT)y + ^ r-n (2) 
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K = py b uf/y S ub- ( 7 ) 



Clearly, formula (HJ is consistent with the log-law (see (@|), and at the same time describes the finite Re modification 
(see (O) taking into account the effect of outer flow geometry (with the parameter m), so it naturally reconciles the 
log-law and power-law debate. All functional forms from to ^ are analytically given, which is a definite step 
beyond the LPR model @. The purpose of this paper is to present analytic derivation of ((T|), employing a novel 
Lie-group formalism. 

The main hypothesis in this derivation is that turbulent fluctuations recover dilation symmetry (discussed below) 
either from the wall or from the center (for channel or pipe) , which is believed to be a general feature of wall-bounded 
turbulent shear flows. Such dilation symmetry manifests in a power law scaling of a relevant length (order) function 
defined using both fluctuation and mean quantities. The introduction of the order function [4] is inspired by the 
observation that a real turbulent flow system (such as a turbulent pipe) differs in a fundamental way from idealized 
system of homogeneous isotropic turbulence or homogeneous shear turbulence by involving a series of statistically 
distinct states of fluctuations. For example, fluctuations near the wall in the so-called viscous sublayer in a pipe 
is very different from that in the bulk flow. In each layer/region, the energy dynamics display different dominant 
mechanisms. For instance, in viscous sublayer, pressure transport balances the viscous dissipation and diffusion, while 
in buffer layer, production (by the mean shear) balances viscous dissipation. As the number of terms in the energy 
equation is finite, the number of balance mechanisms is finite, hence the finite number of layers. 

This multi-state/multi-layer feature is believed to be generic for all systems at a truly far-from-equilibrium state; 
wall turbulence is one of them. Any measure that quantifies the symmetry property of each state and also describes 
transitions between one state to another is thus named order function, in close analogy to 'order parameters' in 
statistical physics [6] , which effectively captures the unknown effects of fluctuation structures on the mean and reveals 
phase transitions (which means from one macroscopic state to another) at macroscopic scales. Our order function 
plays exactly the same role - capturing effects of turbulent fluctuations on the mean and displaying varying scaling 
symmetry in the mean profile. A specific fluctuation field responsible for the vertical momentum transfer is the 
product of streamwise and vertical fluctuations - the Reynolds stress - which is to be described by the mixing length 
(order) function. Our central task here is to characterize the statistical state of each layer using the scaling of the 
mixing length, and then the solution of the mean momentum equation. In the future, a complete deductive theory 
would explain why each layer is necessary and how its quantitative property is derived. 

The paper is organized as follows. In section 2, we establish a formal framework for the Lie-group analysis applied 
to the MME, where an extended phase space including the mixing length and its gradient is introduced. The analysis 
is applied in Section 3 to derive a specific solution of the mixing length for turbulent channel and pipe flows. In 
section 4, we validate the theory using DNS data. Section 5 contains a summary and discussion of future studies. 



II. LIE-GROUP ANALYSIS OF THE INNER AND OUTER MME 



Lie-group analysis is a well-founded method of finding invariant solutions of differential equations |10[ . Recently, 
Oberlack [llj and Lindgren et al use it to derive local scaling laws in the mean velocity, including log-law and 
power-law, for unclosed Reynolds-averaged Navier-Stokes (RANS) equations. When studying multi-point statistical 
symmetry, Oberlack and Andreas [l3| assert that there exist a large number of possibilities for defining invariant 
groups, due to the possible new scaling arising from the fluctuations. On the other hand, empirical data from 
DNS show well-behaved mean profiles, while finding new governing principles to derive it is a challenge. For a 
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symmetry analysis, this amounts to selecting the right invariant groups, which likely correspond to statistically restored 



symmetries by fully developed turbulence, postulated by Frisch 14]. Finding these right symmetries embedding true 
invariant solutions of the RANS equations (i.e. the MME) is an important step forward in the symmetry analysis of 
turbulence. 

A key realization here is that although the RANS equation and its boundary conditions for a channel or pipe 
possess a dilation group (but not a Galilean group), the empirical mean velocity solution seem to break the dilation 
invariance. Even locally, the mean velocity does not have a recognized scaling invariance, except near the wall with 
a rather trivial linear scaling in y + . Specifically, when the RANS equation undergoes a dilation transformation, the 
group invariant of the mean velocity solved by using the Lie-group method is no longer invariant, as can be easily 
checked with any empirical mean velocity profile. Classical Lie-group analysis seemed to meet a bottleneck, yielding 
very few results. This obstacle is overcome here by the introduction of additional variables, namely mixing length 
and its gradient (both order functions), as group invariants. 

Note that in the statistical physics of phase transitions, order parameters are considered to the primary variables 
describing the macroscopic state; we carry out the same tradition: regarding the order functions (i.e. mixing length 
and its spatial gradient) as fundamental variables describing the dilation-symmetry of the underlying fluctuation field 
(i.e. Reynolds stress). In other words, when undergoing dilation transformation, the mixing length and its gradient 
display specific scaling in each layer, and a generalized symmetry relating the mixing length with its gradient describes 
the transitions between the layers. This Lie-group formalism thus enables us to show that there exists a dilation-group 
which leaves the governing equation (i.e. the RANS equation) invariant under the postulated scaling symmetry for 
the mixing length. While the mean velocity is itself no longer invariant under the dilation transformation, the mixing 
length maintains the scaling symmetry. 

This unconventional use of the Lie-group analysis embodies a physical insight about the problem, namely that the 
mixing length and its gradient are the ri ght variables to represent the symmetry property of the fluctuations (e.g. the 
Reynolds stress). As noted by Cantwell [l(|, group theory begins when all the science needed to deduce the solution 
has been established. In wall turbulence, the physics that has been overlooked so far is the existence of the dilation- 
group symmetry displayed by the mixing length in each of layers like sublayer, buffer, bulk and core. A general feature 
uncovered here is that distinct layers do possess symmetries, and there exist appropriate order functions to display 
them. Including order function as additional variables in the symmetry analysis of the RANS equation marks the 
essential difference between this work and all previous symmetry analysis of the RANS equation. 

For a channel/pipe flow, the MME is [15| : 

d 2 U+ _ d(uv) + _ 1 

dy+ 2 dy+ ~ ~R^' ( ' 

where U is the streamwise mean velocity, u and v are streamwise and normal fluctuating velocities, Re T is wall friction 
Reynolds number, and other notions are described above. Integrating along the normal direction yields 

S + + W + = r (9) 

where S+ = dU+/dy+, W+ = -{uv)+ and r = 1 - y+/Re T ). Prandtl's 5] mixing length, see FigfTJ defines a length 
function linking the Reynolds stress and the mean shear stress: 



~~s+~ 

Then, the MME becomes: 

l M 8+2 + S + -r = 0. (11) 
Note that l\ t fully determines the mean shear stress, 



-1 + ^Jah 
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M 



S+(y+;Re T ) = * (12) 



and then the MVP. Prandtl @ assumed a linear, 



m = W + (13) 



which leads to 



s+ « -L « -L (14) 
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FIG. 1: Profile of the mixing length as a function wall distance y + . The solid line indicates the near wall power- law scaling 
with exponent 3/2, dashed line the linear scaling in the overlap region according to 5:]. Inset shows the central power-law 
scaling with exponent —1/2, where £m = t^/Re T , and r = 1 — y is the distance to the center. DNS data from Iwamoto et al. 
[3 at Re r = 650, Hoyas and Jimenez [l7| at Re T = 940, Wu and Moin [3 at Re T = 1142. 



since W + sa 1, and therefore the logarithmic law for MVP 

U+ « - In (y+) + B (15) 

with Karman constant k ps 0.41 and B s» 5.2 15]. The logarithmic law is widely used in engineering models, albeit 
invalid in the near wall region and center flow because £m departs from the linear scaling (|13[) when the wall distance 
y is out of the overlap region; see in Fig[TJ In fact, in the viscous sublayer, Taylor expansion [l5[ yields W + oc y +3 , 
S + ~ 1, thus oc y +3 / 2 . In the core region, according to the central symmetry, the mean velocity profile is parabolic, 
thus S + oc r; which together with W + ~ r, gives £m oc \f\fr. Note that these two power-law scalings are sharply 
different from (|13j) . Although various modifications have been made ever since, such as the damping function by van 
Driest [l9[ and the wall function in the LPR model @ , a derivation of mixing length function with correct asymptotic 
behaviors is not achieved yet. 

The Lie-group analysis presented below is novel in two aspects. First, by substituting in W + = (£^ f S' + ) 2 , we rewrite 
the governing equation ((5J) as: 

* ~ dy+ 2 + Um \ dy+)\ dy+* ) + Mm * m [dy+ ) + Re T ~ ° [ b) 

with wall conditions: U + (0) — ^jy-(O) = 0. Note that ([5]) is an ordinary differential equation, but since we consider 
below a family of solutions for varying Re T and regard Re T as an independent variable, the partial differentiation 
is introduced in (1161) . This is the so-called inner MME (in wall units), and the primary variables are expanded to 
include the mixing length tj, and its spatial gradient = d£\ [ /dy + . Similarly, in central unit r, we have the outer 
MME (£ M = d£ M /dr): 



1 d 2 U+ o2 dU+ d 2 U+ • fdU+\ 
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^=^-^ + 2^^ + 2Wm^J -1 = (17) 

with center condition U + (0) = U+,£m(0) = oo where £/+ is the center line velocity. Note that in (fT6|) and $T7\ . in 
addition to independent variable y + (or r) and Re T , and dependent variable U, we now have an expanded function 
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space including £m and £m- This expansion is essential, since it allows one to formulate specific assumptions about 
the scaling symmetry of the fluctuation fields (or of some specific measure of it, e.g. the Reynolds stress), which then 
allows one to specify the solution manifold. In other words, using £m and £m, we quantify the symmetry behavior of 
the Reynolds stress, and use it to derive a solution for the mean velocity. Our Lie-group analysis not only identifies 
simple scaling symmetries of £m, locally valid in each of layered structures, but also defines a generalized symmetry in 
the prolonged space of £m and £m, describing a universal transition between one layer to another, and then fulfills the 
goal of defining a complete solution for U. The generalized symmetry is a truly hidden symmetry which describes the 
multi-layer structure of wall-sheared turbulence. Since the multi-layer structures are generic to far-from-equilibrium 
systems, our treatment is of general interest for many flows. 



A. Symmetry transformation 



The analysis consists of the following three standard steps. Step one is to define the symmetry transformations. 
Here, the Lie group transformation S e can be expressed as 

t* = $(t, x, e); x* = (*, x, e), (18) 

where the superscript * denotes transformed variables, t — (y + , Re T ) and x — (U + , ) are independent and dependent 
variables, respectively, and $, \& are analytic functions of t, x and a continuous parameter e. In this setting, the 
symmetry transformation satisfies 



t, a?) = O a:*) =0. (19) 

Note that S £ can also be expressed in an infinitesimal form: 

t* = t + ~i(t,x)e+ 0(e 2 ); x* = x + lf(t, x)e + 0(e 2 ), (20) 

where £ = d$/de\ e= o = (£, y +,^Re T ) and if = d^/de\ e= Q = (riu+,i]i+) are the so-called infinitesimals. Then, Taylor 
expansion, ^(t* , x*) = ^(t, x) + eX c (f + 0(e 2 ), yields the following equation for the infinitesimals (llj |: 

^k=°=(t~+t~ +X p ^ = 0, (21) 

where X is the so called Lie-group operator, and X p = ^^+8^+ +7/^+9^+ +Ve+ ®i+ ^ s referred to as the 'prolongation'. 
(|2f p is an over-determined set of linear homogeneous differential equations for infinitesimals, whose general solutions 
can be conveniently given by the software MaplelS, yielding: 

£ y + = aiy + + bi, CRc T =a 2 Re r , 

Vi+ = (-ai + a 2 )£i I /1 1 r lu+ = (2ai - a 2 )U+ + b 2 , (22) 

where coefficients and hi [i = 1,2) denote dilations and translations, respectively. Now, the no-slip wall condition, 
i.e. U + = = at y + =0, leads to b\ = b 2 = 0; thus, near the wall, only the dilations are permitted. Normalized 
by the parameter a\, we have 

e y + = V + , Crc t = (l + 2a)Re T , 

Vt+ = ol£ + m , f] U+ = (1 -2a)U+, (23) 

where a = (—1 + a 2 /a\)/2 is the only remaining parameter specifying the local scaling of the mixing length. In 
addition, the infinitesimals for the derivatives are 77^+ — —2aU + , rjjy + = (—1 — 2a)U + and 77/+ = (a — l)£^j using 
the contact condition. 

Here, we briefly discuss the meaning of the above symmetry transformation. Consider the transformation from the 
state characterized by (y + ,Re T ) to that by (y + * , Re*). The transformed variables following (|2"3"|) are: 

y+* = e c y+, Re* = e (1+2a)e Re T , 

U + * = e^-^ c U+, £\; = e a *£+, £+; = e^H+j. (24) 

We call a symmetry in U + to denote the fact that the actual (empirical) solution of the mean velocity U + at the 
latter state equals the transformed mean velocity, namely U + * . In other words, if U + * does not equal the actual 



FIG. 2: Sketch of the symmetry analysis for a turbulent pipe flow. A real turbulent state A = (y + , Re T , U + , , t^) is 
denoted by a solid line. Consider a dilation transformation (|24[) to a new state A* — (y + * , Re*,U + * , £^ ,£^j ). Compared 
the transformed state A* (denoted by a dot-dashed line) with the real state B = (y + * , Re* T ,U +ES ,£\ ES ,£\ ES ) (denoted by a 
dashed line). Here, the superscript ES denotes empirical solution provided by experimental or DNS data, which is the true 
solution at (y + * , Re*). The symmetry analysis consists in asking: is there any similarity between A* and B? Previous works 
f [ill. H3|) focus on the symmetry in the mean velocity itself, namely demanding U + * = U +ES ; we argue that this symmetry is 
broken, at least in most parts of the flow domain, that is, U + * U +ES . Instead, we postulate the symmetry of the mixing 
length and its gradient, because they represent the restored symmetry [l4| by turbulent fluctuations (e.g. effect of the Reynolds 
stress on the mean shear). 



solution, we say that the symmetry in U + breaks down. Previous works [TT], EH applying the Lie-group analysis to the 
RANS equation were restricted to the case with no symmetry-breaking of U + . But, if one wants to derive solutions 
of the RANS equation from symmetry analysis, the situation of symmetry-breaking needs to be considered, since the 
empirical MVP from either DNS or experimental data do not seem to have a simple scaling form. 

The new idea introduced here consists in recovering the statistical restored scaling symmetry by turbulent fluctua- 
tions in wall-bounded flow with extra variables, namely the mixing length and its gradient, see Fig. [2j By assuming 
scaling symmetry in ttj, we assert that the transformed mixing length, £j^, is the actual solution at the latter state 
corresponding to the independent variables (y + *, Re*). This postulated scaling symmetry leads immediately to a 
power-law solution of £m, which corresponds to the similarity solution of the first kind described below. 

Another similarity solution (of the second kind) corresponds to a case of symmetry- breaking in £m, but a maintained 
scaling symmetry in £m, which is a solution of defect power-law. A third kind of similarity solution is derived under an 
assumption of a generalized symmetry relation between Im and £m, when the scaling symmetries of both quantities 
break down. By including the prolongation of £m, only three possibilities appear; hence, our symmetry consideration 
yields a fairly complete description of the dilation-invariant solutions to the RANS equation. These details are 
presented in section 2.2. 

Step two is to define group- invariants, which are obtained by solving the characteristic equation: dti/d^i — dxi/drji, 
written as: 



dy+ _ dRe T _ d£\ l _ d£\ l 



y+ (l + 2a)Re T a £+ (a-l)£+ 

dU+ _ dU+ _ dt)+ 
~ (1 - 2a)U+ ~ -2aU+ ~ -1 - 2aU+ 

The solutions to (|25[) define following group-invariants: 



(25) 



H = Re T /y +( - 1+2a \ F 1 =£+/y+«, F 2 = ^/y^- 1 ', , , 

Gi = U+/y+^- 2a \ G 2 = U+y +2a , G 3 = U+y+^+ 2a l { ' 

These group-invariants, functions of y+ and Re T , are the so-called similarity variables [Toj . Note that F\ and F% are 
equivalent to the invariant q = y / 1 \fxvJTTZ, defining the similarity variable in a Blasius boundary layer. 
Finally, in step three, we rewrite = in terms of the group-invariants, 

V = G 3 + 2F 2 G 2 G 3 + 2F X F 2 G\ + 1/H = (27) 

Geometrically, ^ = defines an invariant-surface in the variable space of t and x, which is invariant under a dilation 
group of transformation, and the infinitesimals are tangent vectors to the invariant-surface. The so-called invariant 
solution, for example, Q(t,x) = 0, satisfying the invariant surface condition, i.e. A"0|e=o = 0, must be a function of 
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similarity variables. The present work seeks invariant solution in an extended coordinate space including the mixing 
length and its gradient, namely ©(i^.i^) = 0, (e.g. £„+, ri e + and 77;+ ). 

" M "if 

Similarly, the transformed variables for the outer MME are 

r* = e e r, Re* = e (1 / 2+a)e Re T , 

U+* = U+ -e^ 2 - a ^ e (U+ -U + ), t M = e ae i M , t M = e^- 1 >£ M . (28) 

and we derive a set of group-invariants for the outer flow, by performing the above analysis to the outer MME (in 
center units): 

H = Re T r 1 / 2 + a , F 1 =£ M /r a , F 2 =£ M /r a -\ 

Gi = (17+ - U+) /r 3 / 2 -", G 2 = f/+/r 5 / 2 - Q , G 3 = U+/r 7 / 2 - a [ ' 

where • denotes r-gradient; the outer MME in terms of group-invariants is 

Jf = -G 3 /H + 2F 2 G 2 G 3 + 2F 1 F 2 G 2 2 - 1 = (30) 

It is interesting to contrast (|27[) and (130[) with similar forms of the Blasius equation for a laminar boundary layer. 
Following 10], the Blasius equation written in streamfunction ip is 

dip d 2 ip dip d 2 ip d 3 ip ^ 
dy dxdy dx dy 2 dy 3 

It can be demonstrated that this equation is invariant under the following dilation symmetry transformations: 

x* = e e x; y* = e Xe y; iP* = e (1 - x)e iP (32) 
The group parameter A is determined by the free stream condition, i.e. 

dip dip* 

~q \y — 00 — Uoo — COTlst. — — — oo (^3) 



which leads to 



with the following group invariants: 



Substituting them into (13T1) yields 



This equation can be re- written as 



A=l-A; A = 1/2, (34) 



/= r-r T — ; g= , 7= : ( 35 ) 



£1 + 1*1 

dq 3 2 dq 2 



(36) 



Xui + \XX H = (37) 

where Xi represents the ith order derivative of the invariant /. 

The Blasius equation (|3"?| are the equivalent of our equations (|27p and (|30p . but with an important difference: (|3T)l 
is a closed equation, while (|27|) and (|30]) contain two dependent variables (U and £m)- The well-known closure problem 
of turbulence requires the introduction of closure hypothesis. In a symmetry analysis, closure hypothesis comes from 
symmetry assumption. In the previous attempts [ill. [l2|. symmetry is imposed simply on the mean velocity, that is, 
the group invariant of U + is assumed to be constant, i.e. 

Gi(y + ,i?e T ) = Gi(y + *,Re*) = const. (38) 

As shown in ([HI), with this assumption, candidate solutions of power- law and defect power- law can be obtained for 
the mean velocity, 

U+ = G iy +P , or U+ = U+ - Gif*. (39) 
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The power-law solution finds its correspondence with the linear scaling in the viscous sublayer with (3 = 1; the defect 
power law solution, U + = J7+ — G\y® in coordinate y, corresponds to flow solution near the central part of a channel, 
according to Oberlack fT~D . with fitting parameters G\ = 5.83 and (3 = f.69. However, it can be shown that this 
latter solution contradicts the central symmetry condition, J7+ — U + oc r 2 , hence should be regarded as non-physical. 
Therefore, we consider that imposing the invariance on U to derive closure solutions of the RANS equation is not 
successful. 

With the introduction of the new dependent variables, £m and its gradient, the situation is different. First, we are 
given new freedom in constructing the solutions, because the symmetry associated with the fluctuations (included in 
£m) can be treated separately. Secondly and most importantly, we can independently verify if the assumed symmetry 
is true. If Im is not a good variable, we can study a variant of it, until we find the good order function. Thirdly, 
as shown below, we are able to derive rigorously two local solutions near the wall and around the centerlinc of a 
channel/pipe, similarly as deriving the 1/2 exponent in (|34[) for laminar boundary layer. Finally, this program yields 
a general procedure, with the concept of the order function, of performing symmetry analysis for unclosed problems 
involving fluctuations. Below, we show that only three kinds of symmetries are sufficient for determining the full 
profile of mean shear and hence mean velocity according to (|12p . which is highly encouraging. 

B. Three kinds of similarities for mixing length 

When considering the possible symmetries of and t^, three cases naturally arise. 

1. Power-law 

The simplest situation is that the scaling symmetries of tL and are both maintained, as the group parameter 
e continuously changes. This can be sketched by the mapping: 



(y + \ 




(v + * \ 




( = y+* \ 


Re T 




Re* 




= Re* 


{+ 








_ f +*ES 


V + M ) 




{in J 




\ = Z + M ES J 



(40) 



Trans f ormed— variables Empirical — solution 



where the superscript ES denotes the empirical solution. With this symmetry assumption, (|26|) yields constants F\ 
and F 2 : 

Fx = a, and F 2 = ac\. (41) 
Solving (|2"o) . we thus have a power-law solution for t^: 

tii = ay +a (42) 

Note that F2/F1 = a is a key quantity to characterize the solution. Since the scaling symmetry is postulated, the 
actual solution may realize only a specific value of a. Thus, the scaling symmetry assumption provides a candidate 
solution to be validated by empirical data. Specifically, the theory proposes to plot the empirical data in terms of a 
so-called diagnostic function: 

7 - F2/F1 = dln{l+ I )/dln(y + ), or 7 = F 2 /F 1 = dln(£ M )/dln(r). (43) 

If the empirical 7 displays a plateau in a range of y + (or r), then, a local scaling symmetry of is validated, and the 
value of the plateau is thus a which guarantees the dilation invariance of the RANS equation. Such tests carried out 
for DNS channel and pipe flow data indeed show that a = 3/2 in the viscous sublayer, and a = —1/2 in the central 
core. Both exponents can be derived from physical invariants at the boundaries as Q33p for the Blasius equation, see 
Section 3. 



2. Defect Power-law 

The second situation is that the scaling symmetry of breaks down, but that of is maintained. This takes 
place in the bulk flow (for r = O(l)), where the mixing length tends to a finite value of the order of the radius of the 
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pipe (or the half width of the channel). This finite length breaks the scaling symmetry and any dilation invariance. 
Using the outer variable r for the outer MME, we then have the following mapping: 



(v + \ 




(V + * \ 




= y + * \ 


Re T 




Re* 




= Re* T 


£ M 


=> 


0* 
M 


— > 


i u+*ES 













Transformed— variables Empirical — solution 



It can be expressed as 

Fi^ci, and F 2 = c 2 . (45) 

This solution is postulated for describing the bulk flow. Integrating (|29|) once using (|45|). we obtain £m — d +Fir m /m, 
where d is an integration constant and m is an exponent characterizing the power-law of £ m . Furthermore, using the 
boundary condition: lyi — > as r — > 1, then d = —F-xJm. Thus, 

£m = -F 2 (l-r m )/m. (46) 

Note an interesting asymptotic behavior of (|4"6"1) : as r — > 1, Im — > —-^2(1 — r ) = thus we obtain n — —F 2 . In other 
words, the classical Karman constant is a Lie-group invariant of the gradient of the mixing length which is valid not 
only for the overlap region, but for the whole bulk region also. Here, we call (|45l) the similarity of the second kind, 
and PB1 the defect power-law. 

Note that (|46[) is obtained by employing a prolongation in the Lie-group analysis, which deserves a further comment. 
Bluman and Kumci [20] proposed to describe the similarity of the first kind as displaying invariant curve, while the 
similarity of the second kind (with the symmetry-breaking of the original scaling variable) displaying an invariant 
family of curves parameterized by e. In other words, the latter case is characterized by an invariant manifold composed 
of a family of invariant curves. Here, we see a general rule, when an invariance property breaks down (i.e. original 
invariant curve is no longer invariant), another weaker form of the invariance (i.e. a family of invariant curves or 
invariant manifold) would be ensured. This feature can be formulated as a "weak-invariance-principle" which may 
guide the symmetry analysis of complex problems (with substantial symmetry-breaking) . 



3. Composite solution 



Finally, consider the situation when both symmetries in t^j and t^j break down, as sketched below: 

(y + 

Re T 

f + 

V4 J 







( = y+* \ 


Re* T 




= Re* T 


M 




1 f+*ES 


V ft ) 







(47) 



Trans f ormed— variables 



Empirical — solution 



In fact, this happens during a transition from one scaling state to another, when there exist several local scaling states 
in a system composed of a multi-layer structure, with exponents 7/ (/ = 1,2, ...) each characterizing a layer. What 
would happen now? The weak-invariance-principle formulated above applies here, and suggests finding a new form 
which is yet left invariant under the group transformation; this would be then a generalized symmetry which describes 
a continuous transition from one power-law solution, t\l = F(y +l1 to another, t\l* — F^y^ 11 . Let us sketch this 
generalized invariance law using an abstract function G as: 



0(4rVi*) = O 



Trans formed— variables 



e ss 0?+V+*) = o 

Empirical — solution 



(48) 



A specific ansatz is proposed as follows: 



9= ( 7 -7/)-(7//-7/)(^i // M) 



(49) 



where j(y + ) — F 2 /Fi and Fi(y + ) = £~^.f(y + )/y +1 " are variables, while all other symbols are constants. In particular, 
n is a parameter specifying the transition sharpness. (|48[) implies that = is an invariant property under the dilation 
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group transformation; it can be shown that it satisfies the so-called invariant surface condition (see appendix). In 
other words, © = is a generalized symmetry - an invariant surface in the function space of ij^^ttf that is tangent 
to the Lie-group operator, in analogy to the simple power-law and defect power-law solutions. The postulated ansatz 
(|49p imposes a specific invariance form, which would be constructive if (|49j) yields the actual empirical solution. Our 
subsequent studies, reported here and elsewhere, show that it is indeed the case. 
One may wonder how (|4"9")) is motivated? In fact, O = can be rewritten as 

7-7/ / FP \ ™ 

7 11 - 1 1 » (50) 



Hi -71 \ Fx 

For adjacent layers I and //, 7 = F2/F1 as a function of y + , must vary smoothly and monotonically from 7/ to 
7//. As a result, it is easy to verify that (f5TJ)) has the desired asymptotic behavior: 7 goes to 7/ at one end and 

to 7// at the other end. Let y+ = (F{/F* IS \ 1 " 71 ^ denote such a location in the middle of the transition that 
£lj(y+) = llPiyt). Then, for y+ > y+, F x -> F/ 1 , the r.h.s. of {H goes to 1, as expected. For y+ -C y+, 
F 1 /Fi —> {y + /yt) 111 ~ JI 1 an d the r.h.s. of (|50l) goes to zero for p = nfaii — 7/) ^ 1. The latter can always be 
guaranteed by choosing an appropriate n (thus n specifies the sharpness of the transition between two adjacent layers). 
Therefore, (|50|) connects the two asymptotic states of power-law invariance with desired scalings. 
From (|50p . we find an analytical expression for valid across the two layers as: 

e$ I - II) = Fly*<B I -< I ( y + ]y +,p) (51) 

where the transition function B is defined as 

B I -"(y+-yt, P )=(l + (y+/yt) P ) { '" I ~" l)/P - (52) 

It works as a blending procedure to switch between layers, which has been used as an interpolation function in many 
models, e.g. Batchelor [21j used it to describe the transition of the velocity structure function from the dissipative 
to inertial range. We derive it only by a continuity and smoothness argument on the Lie-group invariants. Here, we 
show that a dilation- invariant structure is behind B; it would be intriguing to investigate why this form is so generic 
and what determines the transition sharpness. 



III. THE MULTI-LAYER DESCRIPTION OF TURBULENT CHANNEL AND PIPE FLOWS 

Now, we apply the above analysis to turbulent channel and pipe flows. The notion of the multi-layer comes from 
empirical observations well known in the literature [L5| . For the inner solution, one first meets with the so-called 
viscous sublayer very close to the wall; then, a buffer layer where coherent near-wall structures such as streamwise 
vortices dominate the dynamics [22]. Then is the overlap region between the inner solution and the outer solution. 
For the outer solution, there is a bulk flow, with a core around the center. With the Lie-group formalism developed 
above, these layers can be described analytically and quantified with parameters related to local scaling, transition 
location and transition sharpness. Below, we derive analytical expressions of the mixing length for all of these layers. 



A. Inner flow 



The viscous sublayer is assumed to hold a similarity of the first kind, and is analytically described by a power-law 
scaling of the mixing length as: 

£+° ub = Fl ub y+V 2 (53) 

This result can be justified by a near- wall expansion in the small fluctuation amplitude: u ~ y + and v ~ y +2 such that 
W + ~ y +3 . Since S + s» 1 near the wall, i\ t ~ y +3 ^ 2 , consistent with (|53[) . Note that the symmetry transformation 
itself does not select the exponent 3/2; only further physical consideration related to the boundary constraint is able to 
offer a concrete proposal. Following the argument leading to (p4|) . a similar constraint as (|33|) is S + = 1 at vanishing 
y+, which yields W+/1+ 2 = 1 = W + */(£+;) 2 - (y + *) 3 /(^M ? ■ Thus, a = 3/2 in (23). 

For large y + when (px?)) ceases to be valid, we postulate new layers with the power-law scaling of a distinct exponent, 
and identify them with the buffer-layer and the log-layer. The reason behind this assumption lies in the belief that the 
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dilation-group invariance should be a robust property throughout the whole near-wall region. Indeed, an empirical 
study with the diagnostic function 7 of the DNS channel flow data confirms our assumption with the following results: 

e+w = V 2 (54) 



C 9 = F[° 9 y + = ny+ (55) 

The exponent 2 for the buffer layer is believed to be related to coherent vortical structures, such as streamwise vortices 
whose core is mostly located at y + ~ 25 [22J, which will be studied in details elsewhere. 

Then, the generalized symmetry - ansatz (|50[) - yields a composite formula valid across both the sublayer and buffer 
layer as 

4 (s " 6 - 6u/) =i + M sub B s - b (y + ;vtu b >Pi)> (56) 

where B s - b = (1 + (y+/y+ b ) Pl ) 1/2pi ■ Similarly, 

i+ J b uf-lo g) = e+ Juf Bb - t{y+ . y ^ p2) (57) 

valid across the buffer and log layer, where B b ~ l = (1 + {y + /y bu f) P2 ) _1 / p2 . The test of these formula is presented 
below. 



B. Outer flow 



As mentioned above, the bulk flow is characterized by the finite mixing length which leads to a dilation-symmetry- 
breaking. By assuming that the dilation-symmetry is still maintained in the gradient of the mixing length, we arrive 
at a defect power-law solution, 

i b ^ lk = k(1 - r m )/m. (58) 

Empirical study shows that m = 4 for channel and m = 5 for pipe. We have then a similarity of the second kind. 
Note that ([58]) smoothly match the log-layer = ny independent of m. The bulk solution (|58|) is also a discovery 
of the present theory. A verification by the DNS channel flow data is presented in Fig. 2, where the compensated plot 
of l\i by (|55]l clearly shows a plateau. 

Near the centreline (r —¥ 0), we find an extra layer called core, where the mixing-length diverges as £m oc l/\/r- 
This is consistent with a similarity solution of the first kind: 

gcore = F core r -l/2 ^ 

where the exponent "f core = —1/2 is derived from an invariant central condition 

1 dU+ 1 dU+* , , 

--5— = ~z a „ = const. 60) 
r or r* or* 

which, using J28|), leads to 

3/2-a = 2; a = -1/2. (61) 

This determination is again similar to the argument used to derive (I34[) in the Blasius equation. It is consistent with 
a central symmetry argument: as r — > 0, W + ~ r and S + oc r, hence, £m oc r" 1 / 2 . Furthermore, using (|50l) . the 
Lie-group theory yields a matching from (|59p to a constant mixing length in the bulk flow with scaling exponent zero, 
(-M ~* ^0 = n/m; and the composite solution from the bulk to the central core is thus 

^bulk-core = {fi / m)B b-c {r . ^^3)/^ ( 62 ) 

where B b - C = (1 + (r/r COI . e ) P3 )" 1 / 2p3 and Z c = B b - C (r = 1) is a normalization factor. 
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C. A complete formula using multiplicative rule 

By applying similarity solutions of all three kinds, we have achieved a series of analytical expressions, valid for 
sub-buffer layer (1551) . buffer- log layer (|57p. bulk-core layer. Now we attempt to connect them together to form a 
complete formula for the mixing length profile, i.e.(fT]). 

Following the multiplicative rule, cj) 1 - 111 = (f) 1 - 11 ^ 11 - 111 /<f> Common , is applied to ([56]) and ((57]), yielding a 
matching solution for the inner flow 

£+/" = t+° ub B s - b B b ~ l (63) 

where the factor for the common state cancels. Similarly, applying the multiplicative rule to ([58]) and (1621) . the 

outer solution is 

e Out = e bulk B b-c /Zc (64) 

Finally, applying the multiplication rule to (|63[) and (|64p , we obtain a matching solution for the mixing length profile 
valid in the entire domain: 

f+ _ n+In p+Out I p+Common n+In n+Out l n+log i nr\ 

l M ~ l M M I M M *-M l t M ' V ud J 

which yields (fl}. 

Note that £~^f ut = i?e T £°"* in (|65]) : and p = nytublvtuf m ©• The coefficients in ([I]) can be measured locally and 
independently since they are derived layer by layer: yf ub , y^ u f and r core are the successive transition locations between 
sublayer, buffer layer and core layer, respectively; m — 4 for channel and 5 for pipe represents different bulk; and Pi 
is transition sharpness (i — 1, ...), less sensitive to the prediction of MVP. After all of the parameters are determined, 
([1]) gives a highly accurate mean velocity profile over a wide range of Re (discussed elsewhere). Note that, while 
(TTJ) is established by heuristic argument in |4|, the present work represents a major advance in its derivation, e.g. a 
symmetry analysis of the multi-layer structure in terms of Lie-group language. 

IV. VERIFICATION BY DNS DATA 

We have outlined above a novel Lie-group analysis for determining the mixing length profile. This analysis can be 
subjected to a test by DNS data. In other words, the invariance property of the mixing length function postulated in 
the analysis can go through scrutiny against DNS data which are the true Navier-Stokes solutions. Currently, there 
exist several DNS datasets for channel and pipe flows for Re T up to thousands. We have selected two channel flows 
from Iwamoto et al. [lj| at Re T = 650 and Hoyas and Jimenez [ljj at Re T = 940, and one pipe flow of Wu and Moin 
[l8j at Re r = 1142, which are believed to be well-resolved. The data confirm the multi-layer similarity law for the 
mixing length, as below. 

Fig. [3] shows a clear evidence for the newly discovered bulk flow structure, 1 — r rn (m = 4 for channel and 5 for 
pipe), using a compensated plot of the observed mixing length (e.g. divided by 1 — r m ), in which each of the inner 
layers and the core (divergent) layer are also very visible. The existence of the central core is due to a switch of the 
kinetic energy generation mechanism, from mean shear (production) to turbulent transport. Close to the centreline, 
the mean shear approaches zero, and finite energy dissipation is balanced by turbulent transport. To this end, a 
clear four-layer picture for turbulent channel and pipe flows is fully confirmed! The thickness of each layer is an 
important parameter, sketched in Fig. [3] (the demarcation for the core layer is sketched by y^ ore = (1 — r core )Re T , 
hence increasing with Re), whose accurate determination from DNS data will be discussed elsewhere. 

V. CONCLUSION AND DISCUSSIONS 

In this paper, we demonstrate that a recently proposed multi-layer formula (JIJ for the mixing length [4| can be 
derived by a Lie-group symmetry analysis for the inner and outer MMEs, expressed not only in terms of the mean 
velocity, but also the mixing length and its gradient - the latter representing effects of fluctuations (Reynolds stress). 
Under the assumption about the symmetry of the mixing length and its gradient, three kinds of similarity solutions 
are derived. DNS data soundly validates our description, hence postulates. 

A feature of the new theory is that it expresses the closure solutions in terms of three sets of parameters: scaling 
exponent, transition location (thickness in layered structure) and transition sharpness. These parameters are unique 
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A 




FIG. 3: The mixing length function from DNS data divided by the bulk flow solution reveals a four-layer structure, i.e. viscous 
sublayer, buffer layer, bulk zone and core layer, separated by yf ub , y buf and r COTe . 



to the dilation-group invariance theory, which are, in general, functions of Re, but at large Re, they are believed 
to be Re- independent. Regarding the i?e-similarity, note that in previous works [Tl], [T^j , the Re (viscosity) effect is 
typically neglected by assuming a i?e-independent solution of the mean velocity, which is obviously invalid because of 
the complex matching between the inner and outer solutions. This problem is solved here by writing down separately 
the inner and outer MME, and the i?e-similarity is guaranteed by involving Re explicitly in the transformation. Aa 
large enough Re, the i?e-similarity may be translated to .Re-independence; in this case, the Lie-group invariant F\ 
becomes truly invariant: 

F 1 (y+,Re T ) = F 1 {y +m ,Re*) = Fx{y + *, Re T ) = const. (66) 

where the first equation stems from the symmetry transformation, and the second from i?e-indepcndcncc, as mentioned 
before. Empirical study in a turbulent pipe shows that the Re independence only holds for Re T > 5000, while a slow 
increase of the parameters y^ uf and r core is found below this critical Re. The details will be discussed elsewhere. 
With a model of i?e-dependence of the parameters, the theory ([1]) predicts the mean velocity profile for all Re. 

What does the symmetry analysis using the Lie-group formalism add to our understanding more than making a 
direct postulate of power-law or other ansatz for £m? The answer is that a Lie-group formalism guarantees that the 
RANS equation remains invariant under the (dilation) group of transformation, hence any ansatz used here yields a 
similarity solution. In other words, the family of solutions based on power-law, defect power-law and the generalized 
invariance ansatz of £m form a set of invariant solutions of the RANS equation, when £m is treated as one dependent 
variables. Future study would clarify why £m is so special to possess this invariance; once it is done, a new avenue of 
solving the closure problem of turbulence would be opened. 

Thus, we speculate to have formulated a general framework for turbulence closure problem: searching for group- 
invariant properties of relevant length order functions representative of the fluctuation structures (e.g. the mixing 
length and some other length order functions), and then solving for the mean velocity (which has a symmetry-breaking) 
from the MME. This procedure can be extended to other mean quantities by introducing relevant order functions 
which possess the required symmetry and guarantee the invariance of the relevant balance equations. Ultimately, 
whether such closure solutions are true solutions need to be verified against empirical data. If verified, the length 
order function (and its gradient) is indeed the right quantity, since the order function is, by definition, the quantity 
displaying distinct symmetry. Hence, we achieve a logically consistent framework. 

Note that the Lie-group formalism is a necessary and powerful tool for the present success. Only with the present 
formalism, an analytical expression for the mixing length is obtained, for the first time, wholly motivated from 
symmetry consideration, i.e. dilation-group invariance. While power-law form for the mixing length can be motivated 
from simple scaling arguments, the present formalism accomplishes two additional aspects which would be difficult 
otherwise: a defect-power law form for the bulk flow and a series of matched forms derived from the generalized 
invariance ansatz. Note that the ansatz requires only simple assumptions about the continuity and smoothness for 
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describing the variation of the local group invariants. Finally, it is remarkable that there exist only three distinct 
classes of similarity solutions, associated with the dilation group of the mixing length with its prolongation. 

Finally, we note that the mixing length is only one example of the order functions, specifically representing the 
effect of the Reynolds stress in the mean momentum transfer. The complexity of turbulent flows is that they generally 
encompass a hierarchy set of order functions, each displaying a specific type of turbulent transfers. One will find other 
order functions responsible for turbulent transfer of kinetic energy, temperature, etc. O ur p reliminary study has 
already identified several other order functions in compressible turbulent boundary layers [24j, rough pipe [251 ] ■ and 
Reyleigh-Benard convection. Thus, the concept of the order function is general, applicable to many other flow systems, 
and the Lie-group formalism developed in the present work is also of general interest. The present analysis opens a 
promising direction for modeling the mean profiles in a wide class of turbulent wall-bounded flows. 
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Appendix A: Discussion on the generalized symmetry (|49[) 

Note that (|49p can be written as a function of similarity variables: Q(Fi, F2) = F2/F1 — d\ — d%Fi = 0, where d\ 
and d% are both constants. Substituting F\ — tt I /y +a , F2 — ^Xr/v a into it, yields: 

Q(F U F 2 ) = V + t + M /t M -dx- d 2 (ii I /y +a r. (Al) 

As the Lie-group operator in the variable space (y + , , ^m) 1S: 



X = y+ — +a£+— r + (a-l)^+_ rr . (A2) 
y M 0*+ V 1 M d£+ I V ; 



Hence, 



XQ = X(y+£+/e+ -d 1 - d 2 (£+/y+«r) 

= ~ ay+ih/tb + (a- l)y + l + M /l + M ) - d 2 (na(l+ - na(4/ffl (A3) 

= 0. 



